Dilepton yields from Brown-Rho scaled vector mesons including memory effects 
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Brown-Rho scaling, which has been strongly discussed after recent NA60 data was presented, is 
investigated within a nonequilibrium field theoretical description that includes quantum mechanical 
memory. Dimuon yields are calculated by application of a model for the fireball, and strong modifi- 
cations arc found in the comparison to quasi-equilibrium calculations, which assume instantaneous 
adjustment of all meson properties to the surrounding medium. In addition, we show results for the 
situation of very broad excitations. 



o : 
o . 

(N ■ 

wo: 
s ■ 

< 



> 

(N 

o 

oo 
O 

o 
i 

a: 
43 



X 
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Motivated by the recent discussion on whether Brown- 



Rho scaling [1| is 'ruled out' by the NA60 data, pre- 
sented at Quark Matter 2005 Q, L| , we calculate dilep- 
ton yields from p-mesons with dropping masses within 
the nonequilibrium quantum field formulation introduced 
in p|. First quantitative but simplified calculations 
within this work have shown that memory effects, which 
are neglected in quasi-equilibrium calculations |5j and 
0,0, U0|, have an important influence on the final dilep- 
ton yields, especially for dropping mass scenarios. This 
lead to the conclusion that an analysis of which scenario 
describes the data correctly demands the proper inclusion 
of memory effects. The aim of this work is to clarify how 
the inclusion of memory affects the dimuon yields from 
Brown-Rho scaled vector mesons. The validity of stan- 
dard equilibrium calculations depends on the strength of 
these effects. 

We consider two mass parameterizations: The temper- 
ature and density dependent one used by Rapp et al. |Hj 



i* = m p (l - 0.15n B /n ) [l - (T/T c ) : 



,0.3 



(1) 



assuming a constant gauge coupling g in the vector meson 
dominance (VMD) -coupling, and one motivated by the 
renewed version of Brown-Rho scaling discussed in llll 
0: 



in, 



m p (l - 0.15 n B /n ), 



(2) 



with a modified gauge coupling g* , in such a way that g* 
is constant up to normal nuclear density tiq , while from 



then on m*/g* is taken to be constant 13]. This parame- 
terization is not valid close to T c because it does not take 
the meson mass to zero at the critical point. However, 
Brown and Rho argue 0, E3 that because lattice cal- 
culations show that the pole mass of the vector mesons 
does not change appreciably up to T = 125 MeV, the pa- 
rameterization using [1 — (T/T c ) n ] d (with positive d and 
integer n) overestimates the mass shift. They find tem- 
perature dependent effects to be an order of magnitude 
smaller than the density dependent effects, and hence 
suggest to concentrate on the density dependent part. 



We also wish to study the effect of finite memory for a 
mass shift following this parameterization and ignore its 
shortcomings for the moment. Furthermore, Brown and 
Rho point out that due to the violation of VMD, which 
accounts for most of the shape of the dilepton spectrum 
(see [HEII and 0). the 

overall dilepton production in 
dense matter should be reduced by a factor of 4 com- 
pared to Rapp's calculations, which we do not take into 
account here. 

We compare the usual approach [J^, where the dilep- 
ton rate is given by the well known equilibrium formula 



dN 
d 4 xd 4 k 



(r,fc) 



2e 4 



(2tt) 5 g 2 



■£(M) 



x (T(t), ko) 7tA p (t, k), 



(3) 



with invariant mass M 2 = 
6(M 2 



1 



k 2 = k 2 - k 2 and C(M) = 
— 4m 2 ), to the nonequilib- 



2m, \ 4m? 

rium formalism, in which the propagators of the p-meson 
and the virtual photon are calculated using the general 
nonequilibrium formulas. The dilepton rate for a spa- 
tially homogeneous but time dependent system is given 
by Q: 



dN 
d 4 xd 4 k 



(r,k) 



2e 2 
: (2tt) 5 

x 



M 2 C(M) 



to 



diiD< T (k,r,¥)e ik ^ T -^ 



(4) 



In its derivation we treated transverse and longitudinal 
modes equally, which is adequate for our purposes. D< T 
is the transverse virtual photon propagator, and satisfies 
the generalized fluctuation dissipation relation 



IT 



ret < /-)adv 



D rct q < D 



(5) 



with all time variables and integrations, and further in- 
dices implicit. The relation follows directly from the 
Kadanoff-Baym equations when terms corresponding to 
initial conditions at time to are omitted. This can be done 
when the system is given enough time to reach its initial 
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state before we start the actual calculation of the yield 
0. In this case the £) rct / adv are the free retarded and 
advanced virtual photon propagators while q < — II < , the 
photon self energy, follows VMD by 

„m* 2 (ti) „ rn* n 2 {t 2 ) 
^(h,t 2 ) = e 2 -^D< T ( tl ,t 2 )-^. (6) 

Note that the couplings are to be taken at the correct ver- 
tices, separated in time. The free photon propagator has 
to be exponentially damped in order to reduce large con- 
tributions from early times, which make the extraction 
of the higher frequency structure numerically problem- 
atic 0. 

The meson propagator D< T also satisfies (J5J, with 
q < = S < , the meson self energy, and the retarded and 
advanced meson propagators, which obey the equation 
of motion that in a spatially homogeneous and isotropic 
medium is given by 

- f 1 di"S^(k, t u t)D p %(k, i, t 2 ) = 6(h - i 2 ), 
Jt 2 

(7) 

and Dfj(k,t u t 2 ) = D^(k,t 2 ,*i)- For calculating £< 
from S ret we assume local equilibrium and introduce the 
fireball temperature to the calculation. The integral over 
past times in Eq. as well as the time integrations in 
Eq. © have encoded the finite memory of the system. 
In we show that in order to describe a heavy ion re- 
action it is essential to retain this dynamic information. 
This is due to the time scales for the adaption of the me- 
son's spectral properties to the evolving medium being of 
the same order as the lifetime of the regarded hadronic 
system. In this situation the evolution is non-Markovian 
and quantum mechanical interferences among past time 
contributions become important |J| . Therefore a gradient 
approximation can not properly describe the situation in 
a heavy ion collision and remarkable differences between 
the dynamic and the Markovian yields occur. 

Following 0, H, E|, we use a fireball model with a 
cylindrical volume expansion in the ±z direction: 

V(t) = {z + v z r)ir{r + 0.5 a ± r 2 ) 2 , (8) 

where zq is equivalent to a formation time, v z = c is the 
primordial longitudinal motion, ro = 5.15 fm is the initial 
nuclear overlap radius, and a± = 0.08c 2 /fm is the radial 
acceleration. We start the calculation at T = T c = 175 
MeV (at r ~ 3 fm/c) and set the meson self energy 
S < to zero for temperatures below chemical freezeout at 
T = 120 MeV, thus ending the production of p-mesons. 
The hadronic phase, for which we calculate the dimuon 
production, lives for about 6 fm/c. The time dependent 
temperature and density, determined via the fireball evo- 
lution, yield the time dependence of the p-meson mass 



following parameterizations and J5J). We include ad- 
ditional broadening b y a bout 100 MeV as it was done 
in Rapp's calculation |lOl| shown in In the dynamic 
case the p-meson propagator and the photon propaga- 
tor follow Eq. (JSJ after a Fourier transformation of the 
self energy into the two time representation. The rate 
is then calculated using Eq. (0} , whereas in the quasi- 
cquilibrium case the rate at a certain time follows di- 
rectly from Eq. using the spectral function and Bose- 
distribution corresponding to the present mass and tem- 
perature. The calculation is done for each momentum 
mode between k = and k — 1.5 GeV and integrated 
over momenta. We show the comparison of the dynamic 
and the quasi-equilibrium calculation for the two param- 
eterizations in Figs. ^andE! We find a strong enhance- 
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FIG. 1: Dimuon yield from decaying p-mesons calculated out 
of equilibrium and within the usual static approximation us- 
ing parameterization Q . A difference of about a factor of 4 
is visible in the regime between 400 and 800 MeV. 
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FIG. 2: Dimuon yield from decaying p-mesons calculated 
out of equilibrium and within the usual static approximation 
using parameterization The most significant difference 

(about a factor of 2-3) is visible between 600 and 800 MeV. 

ment as well as modification of shape in both cases. For 
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the density and temperature dependent parameterization 
the dynamic calculation yields about a factor of 4 more 
dileptons between invariant masses of 400 and 800 MeV 
as compared to the Markovian case, while the other pa- 
rameterization shows an enhancement of about 2-3 times 
between 600 and 800 MeV in the dynamic calculation. 

The reasons for the strong differences are the follow- 
ing: First, the meson mass effectively approaches its 
vacuum value more slowly in the dynamic calculation. 
Lower masses are enhanced due to the Bose factor and if 
the spectral function has a lower mass for a longer time 
enhancement is greater. Second, the memory of higher 
temperatures increases this enhancement for all masses. 
Finally, the modified coupling in Eq. © suppresses the 
early contributions. This means that for the static case 
all low masses in the spectral function and high tempera- 
tures are suppressed. However, due to the finite memory 
in the dynamic case, these low masses and high temper- 
atures still contribute at later times when the coupling 
is larger, meaning that they are less suppressed by the 
VMD-coupling. 

In the first case, besides the overall increase in yield, 
we find a significant enhancement in the dynamic calcu- 
lation for invariant masses around the p-vacuum mass. 
Comparing to NA60 data (shown vs. the equilibrium 
Brown- Rho calculation in e.g. the enhancement 

would improve the situation for the Brown-Rho scenario, 
which in the equilibrium calculation strongly underesti- 
mates the data in this regime. The second parameteriza- 
tion leads to even more weight around the vacuum peak 
in the dynamic case and would also improve the agree- 
ment with the data. 

Having investigated pure dropping mass scenarios, we 
now combine Brown- Rho-scaling using parameterization 
with strong broadening and coupling to the excitation 
TV* (1520). This is achieved by using the self energy 



ImE ret (r,a;,k) 



x 9i 



P(t) ( fRNp 

3 V m p 
L> 3 Er R { T ) 



{uJ 2 _ Lj^Z. _ £2)2 + (r fl ( T ) W )2 



-U)T(l 



(9) 



with E = y m R + k — mjv and wir and tun the masses 
of the resonance and nucleon, respectively [lij . Tr is the 
width of the resonance, which in vacuum is 120 MeV, 
and gi the isospin factor. Both widths T and Tr in- 
crease substantially in the medium (r — > 700 MeV and 
r_R -» 400 MeV) 0, E3 - Tnc transverse and longitudinal 
components of the self energy © have been combined, 
following our approximation of treating both contribu- 
tions equally. Note that we are always using the vacuum 
p-mass in Eq. to prevent the coupling to the reso- 
nance from becoming infinite. Typical resulting spectral 
functions are shown in Fig. [3] 
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FIG. 3: Medium modified spectral function using self energy 
© with (Medium 1: m = 300 MeV) and without (Medium 
2) mass shift at normal nuclear density and widths of F = 700 
MeV and r_R — 400 MeV. The p-vacuum spectral function is 
shown for comparison. 



The result of the calculation combining Brown-Rho 
scaling, broadening and coupling to the iV*(1520) is pre- 
sented in Fig. 2| Here, the differences between the 
Markovian and the dynamic calculations seem less signif- 
icant as compared to the mass shift scenarios with only 
marginal broadening (cf. Figs. ^and^. This can be ex- 
pected since the memory time scales approximately with 
the inverse width in the system (see |4J). The very broad 
spectral function adjusts faster to the medium changes. 
However, around the p-vacuum mass the dynamic calcu- 
lation still leads to around three times more yield than 
the Markovian one. 
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FIG. 4: Dimuon yield from decaying p-mesons calculated out 
of equilibrium and within the Markovian approximation using 
parameterization Q for the meson mass and coupling to the 
iV*(1520) resonance by means of the self energy 



For comparison we present the calculation for a sce- 
nario with coupling to the resonance and broadening but 
without dropping mass (Fig. |SJ). The major difference to 
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the scenario with dropping mass is the overall increase 
of the yield. This is due to the fact that in the dropping 
mass case, spectral weight is shifted below the threshold 
of twice the muon mass at early times. For the case with- 
out mass shift the difference between the dynamic and 
Markovian calculations is less than 30% for an invari- 
ant mass up to 650 MeV (above that differences reach 
about 50%) This is again due to the very broad spec- 
tral function, for that memory effects become less signif- 
icant. Furthermore, the time scale for adjustment to the 
medium Q also depends on the pole mass of the spectral 
function. Lower masses mean larger memory times, such 
that for dropping mass scenarios the effects are naturally 
larger. 
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FIG. 5: Dimuon yield from decaying p-mesons calculated out 
of equilibrium and within the Markovian approximation using 
coupling to the A r *(1520) and broadening of both the vacuum 
and the resonance peak to 400 MeV. Maximal differences be- 
tween the two calculations are about 30%. 



In summary, we have shown that the inclusion of fi- 
nite memory by calculation of the dilepton yield within a 
noncquilibrium field theoretical setup, using the realtime 
formalism, yields moderate to strong differences when 
comparing to dilepton yields, calculated assuming instan- 
taneous adaption of all meson properties to the medium. 
Particularly, we have compared the results for dropping 
mass scenarios for dimuon yields in the NA60 scenario, 
which is modelled by an expanding hreball. An enhance- 
ment by about a factor of 3-4 around the p- vacuum mass 
and significant differences in shape have been found for 
the case of dropping mass with moderate broadening by 
100 MeV. Inclusion of coupling to the iV*(1520) reso- 
nance and significant broadening strongly modifies the 
shape of the calculated yields, especially in the dynamic 
calculation. However, the yield is still enhanced and 
more weight lies around the p-vacuum mass in the dy- 
namic calculation. This is important for the comparison 
to the NA60 data, where the Markovian calculation us- 



ing Brown-Rho scaling lead to too few dimuons in that 
region |3J, . A scenario without mass shifts but includ- 
ing strong broadening Q is less affected by the inclusion 
of memory. The dynamic calculation differs from the 
Markovian one by maximally 50%, and less than 30% in 
most of the regarded range. 

We conclude that memory effects must not be ne- 
glected in precision calculations of dilepton yields from 
relativistic heavy ion collisions. Every dropping mass 
scenario we have investigated shows a significant depen- 
dence on whether instantaneous adaption to the medium 
or the full quantum mechanical evolution are considered. 
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We will refer to this as the Markovian approach, because 
memory is completely neglected. 

Explicitly, we use the replacement D t ^(t — ti) = (r — 
ti) — > (r — fi)e _A ' T_tl ' and an analogous one for the 
advanced propagator. Then we renormalize the result by 
multiplication with (lj 2 + A 2 ) 2 /lj 4 (see 3 for more de- 
tails). Note that a variation of A between 200 and 800 
MeV leads to a change in the resulting yield of up to 5% 
in the interesting region between 0.4 and 0.8 GeV. Dif- 
ferences very close to twice the muon mass reach about 
15%. 



